A PURE GEOMETRIC MODEL FOR A MAIN SEQUENCE STAR 
M.I.Wanas*'^ & Samah.A.Ammar*'^ 
Abstract 

A pure geometric field theory is used to construct a model for a spherically sym- 
metric configuration. The model obtained can be considered as a pure geometric 
one in the sense that the tensor describing the material distributions is not a phe- 
nomenological object, but a part of the geometric structure used. A general equation 
Q> ' of state is obtained from, and not imposed on, the model. The solution obtained 

■ shows that there are different zones characterizing the configuration: a central ra- 

i diation dominant zone, a convection zone and a corona. The model can represent a 

^ I main sequence star. The present work shows that Einstein's geometrization scheme 

^ ' can be extended to gain more physical information within material distribution. 

1—5 

in 

^ : 1 Introduction 



Gravity plays an important role in the structure and evolution of macroscopic objects. It 
affects many physical parameters of such objects, e.g. density, pressure and temperature. 
bX), In relativistic theories of gravity, especially general relativity (GR), the relation between 

gravity and such parameters are given by the field equations of the theory. In the case of 
^ _ GR, the field equations within the material distribution are given by [1] 

00 ; 1 

O '. R^lu - T^g^ipR = -i^Tf,„. (1.1) 

00 ■ ^ 

'vh I where T^i, is a tensor giving the physical properties of the material distribution, while 

^ I the L.H.S. of this equation gives the geometric description of the gravitational field. In 

orthodox GR, written in the context of Riemannian geometry, the L.H.S. of (1.1) is a 
O ' pure geometric object while the R.H.S. is a phenomenological one. 

>■ . It is well known that when using (1.1) to study problems concerning structure and 

^ i evolution of large scale objects, one needs to impose an extra condition (an equation of 

5-1 I state) in order to solve the problem. This is done because the field equations (1.1) are, in 

general, not sufficient to get the unknown functions of the problem. This scheme, of im- 
posing an equation of state, leads some times to unsatisfactory results (e.g. Schwartzschild 
interior solution [2]). The solution, in this case, can be considered as a mathematical one, 
without any physical applications. Reasons of this problem may be caused by the R.H.S. 
of (1.1), which is not a part of the geometric structure used, the Riemannian geometry. 
It is preferable, if we use a theory in which the material-energy tensor is a part of the 
geometric structure used. Since Riemannian structure is just sufficient to describe gravi- 
tational fields, one needs a more wider geometric structure, than the Riemannian one, in 
order to describe the material distribution. 
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The aim of the present work is to use a pure geometric field theory, i.e. a theory in 
which the tensor describing the material-energy distribution is a part of the geometric 
structure, in applications. For this reason we give, in section 2, a brief review of a 
geometric structure, with simultaneously non-vanishing curvature and torsion. In section 
3 we review briefiy the field equations of a pure geometric theory written in the context 
of the structure given in section 2. In section 4 we solve these field equations in the case 
of spherical symmetry, restricting it to a special case of a perfect fluid. In section 5, we 
discuss the model obtained. Some concluding remarks are given in section 6. 



2 The Underlying Geometry 

In the present section we are going to review, briefiy, a geometric structure, more wider 
than the Riemannian one. This structure is a version of "Absolute Parallelism" (AP)- 
geometry ( cf.[3]), in which both curvature and torsion are simultaneously non- vanishing 
objects [4]. The structure of the conventional AP-geometry is defined completely (in 4- 
dimensions) by a tetrad vector field such that, its determinant A* = || A^|| is non- 

i i 

vanishing and, 

\^\u = K- (2.1) 

i i 

where A^ are the contravariant components of A^. Einstein's summation convention is 

i i 

carried over repeated indices wherever they exist. Using these vectors, one can define the 
second order symmetric tensor, 

9^lu ~ A/j Xui (2-2) 

i i 

which is clearly non-degenerate. This symmetric tensor can be used to play the role of 
the metric of Riemannian space, associated with the AP-space, when needed. 

Connections, Curvatures and Torsion: 

The AP-space admits a non-symmetric affine connection r"^,^, which arises as a 
consequence of the AP-condition (cf. [3]), i.e. 

Am = A/i,i/ — P" Aa = 0, (2.3) 

where the stroke and the (-I-) sign denote tensor differentiation as defined 0. Equation 
(2.3) can be solved to give, 

r%. = A"A;.,.. (2.4) 

It 

Since P"^,^ is non-symmetric as clear from (2.4), then one can define its dual connection 

[3] as, 

f = r%,. (2.5) 

Also, Christoffel symbol could be defined using the metric tensor (2.2). So, in 
the AP-geometry one could define four linear connections, at least. The non-symmetric 



^ we are going to use (,) for ordinary partial differentiation and (;) to represent covariant differentiation 
using Christoffel symbol 
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connection (2.4), its dual (2.5), Christoffel symbol defined using (2.2), and the symmetric 
part of (2.4). 

The curvature tensors, corresponding to the above mentioned connections, respec- 
tively, are [4] 

T\4-a <^ -pa -pa i pt pa pe -pa (O 

. iJ,i>a . fj,a,i> . iJ,v,(j ' . fia . eu . jxv . ecri \ ) 

-\jfa 4^ "pcK -pa I -pe -pa -pe -pa (n rv\ 

. nua . iJ,a,i> . iJ,v,cr ' . ^la . ev . ij,v . eat \ ) 

ixva — . {ti(T),v ~ ^ . (/ii/),(T + i . (/i(T)^ . {ev) ~ ^ . . (etr) ■ l^-^J 

Due to the AP-condition (2.3) the curvature tensor given by (2.6) vanishes identically, 
while those given by (2.7), (2.8) and (2.9) do not vanish. 

Using the non-symmetric connection one can define a third order skew-tensor, 

A a 4^ "pcK -pa -pa -pa a a /o i r\\ 

. ^IV . //I/ . UfJ, . U/J, . /XU . \ 

This tensor is the torsion tensor of AP-space. One can define another third order tensor 
viz, 

= A" A;.;.. (2.11) 

i i 

This tensor is called the contortion of the space, using which, it can be shown that 

rv = 7.V + C}' (2-12) 

where j^ua is skew- symmetric in its first two indices. A basic vector could be obtained by 
contraction, using any one of the above third order tensors, 

C, = A\a = 7.V- (2-13) 
Using the contortion, we can define the following symmetric third order tensor, 

= + (2-14) 

The skew-symmetric and symmetric parts of the tensor j'^^^ are, respectively 

7.°[H = r"[^.] = Ia\^, (2.15) 

7."(,.) = (2.16) 

where the brackets [ ] and the parenthesis ( ) are used for anti-symmetrization and sym- 
metrization, of tensors, with respect to the enclosed indices, respectively. 
Tensor Derivatives 

Using the connections mentioned above, one can define the following derivatives [4] , 



A'+i."^' A% + '^'aA", (2.17) 
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def 



AH def j^fx Aa 



(2.18) 

(2.19) 
(2.20) 

The following table is extracted from [3] and contains second order tensors that are 
used in most applications. 

Table 1: Second Order World Tensors [3] 



Skew- Symmetric Tensors 



Symmetric Tensors 



e 45.f a 

S/ii^ — Jul/. 



def A Q, 

def 



- C V — C M 



def 



.i^e /.a/i 



def 



d f ' ~^ 

^av — d + C p. 

— -r dsf CK € I CK 6 



de_f £ Q, 
def £ Q, 



a 



def 



/if 



where A" „ = A"*"^ ^ . It can be easily shown that there exists an identity between 2nd 
order skew-tensors, which can be written in the form (cf.[3]), 

Vp^v + - Xixv = 0. (2.21) 
A useful relation between the torsion and the contortion is given by [5] 

IffMau '^(^■^/xau Aj/y^Q, A^y^^). (2.22) 

We see from the Table 1 that the torsion tensor plays an important role in the structure of 
AP-space. All tensors in Table 1 vanish if the torsion vanishes (see (2.22)). The structure 
that is used in the present work is characterized by the dual linear connection (2.5), its 
corresponding torsion (2.10) and curvature (2.7) which are, in general, simultaneously non- 
vanishing objects. So, it is clear that this structure is more wider than the Riemannian 
one and is of the Riemannian- Cartan type. 
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3 The Field Theory Used 



In this section we are going to review briefly the field equations of a pure geometric field 
theory constructed by the authors [6], following a procedure similar to that used in con- 
structing the field equations of GR, in the context of the geometric structure given in the 
previous section. Since this structure is more wider than the Riemannian one, as shown in 
section 2, it will be shown that the tensor representing the material distribution is a part 
of the geometric structure, as expected. The method used to derive the field equations is 
that of Dolan and McCrea [7] variational method His used . 

Field Equations 

The set of field equations, to be used, has the form [6] 

= 0. (3.1) 
This set can be written explicitly in the form. 



+7."%7.\. - r.^'l^. - 2C„7."^ = 0, (3.2) 



where = g^^'N^^. 

To discuss the physical consequences of the geometric set (3.2), it is convenient to 
write it in the covariant form as, 

= -2C + Ng,„ - + 27^ , + 27^ ,^7'^,, 

+ 1: uul'.e - '^C.l\. = 0, (3.3) 



where Sy„ = g^^,S^^. 



The Symmetric Part of 8,^^ 

The symmetric part of Sua is defined as usual by. 

Substituting from (3.3) into the above definition and using the symmetric tensors of Table 
(1), we can write 

S{^ua) = ~2GucT - gua^ + + 2 UJua " <t>u<T = 0, (3.4) 

which can be written in the, more convenient, form 



^ Since Dolan-McCrea method is not published, so the reader may refer to [8] for more details about 
its use in the AP-geometry. 
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where, 

Tta =^ l^i^ya - l^^va + ^va " -j^gva^^- (3.6) 

From (3.5) it is clear that, 

T*''% = 0. (3.7) 

This imphes conservation (since the vectorial divergence of the left hand side of (3.5) van- 
ishes identically). Then the tensor (3.6) can be used to represent the material distribution 
in the theory. 

The Skew- Symmetric Part of Spa 



The skew part of the field equations (3.5) is given by 

which can be written, using the skew tensors of Table (1), as 

= Xva - Vi^a = 0. (3.8) 

Now using the identity (2.21), the skew part of the field equations can be written in the 
form, 

e.a = 0, (3.9) 



4 Solution With Spherical Symmetry 

The tetrad vector field, which defines the structure of an AP-space with spherical sym- 
metry, defined by using the coordinate (a;° = t,x^ = r, = 9,x^ = ip) , can be written 
in the form [9] 







V 



Dr 
B sin 9 cos cp 

sin sin 9? 
BcosO 





y cos 9 cos if 
^ cos 9 sin (fi 
—sm9 



\ 

-B sin (p 



B cos if 
r sin 6 





(4.1) 



where A, B and D are function of r only. We are going to take D = which implies 
the vanishing of all skew tensors that appear in the structure of the field equations (3.1). 
This guarantees the absence of interactions, other than gravity, if any. In this case (4.1) 
will reduce to 



( ^ 




V 





B sin 9 cos ip 
B sin 9 sin (p 
Bcos9 







- cos( 

r 
B 



cos p 
^ COS 9 sin (fi 
^sm9 



Consequently, using (2.1), we get 













~B sin If 

r sin 6 
B cos ^ 

rsinO 









(4.2) 



Am 



sin 6' cos p cos 9 cos p sin ^ sin p 
^ sin 9 sin ip cos ^ sin <^ sin 9 cos (p 
V ^cos9 ^sm9 



(4.3) 
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Using definition (2.2) and the tetrad (4.3) wc get, 



and then, 



f 


1 

A2 


















1 


















r2 


















r2 sin2 6 




S2 



(4.4) 














\ 





















B2 







V 








52 

r2 sin2 6 


1 



(4.5) 



From above tetrad and metric we can evaluate second order tensors of Table 1, necessary 
for solving the field equations (3.3), in the present case. 



Second Order Symmetric Tensors 

The symmetric tensors necessary for the field equations (3.3), using the definitions of 
the tensors given in Table 1, have the following non- vanishing components. 



00 



^11 = -4 



BB'A' 



d2A" I 02^- o 



?2 A'2 



A3 



22 



Br 

r'^B B 

' B ' B 



, "033 = V'22 sin^ 9 



and. 

Also, 
and. 



(P22 = 52 



-4 



A^B 
S'2r2 



9A'2 



- 2 



A'B'r^ 



~AB^ 



, ¥'33 = ^^22 sin^ ^ 



R 
R 

def 



00 — 


B'A'B 
A3 


B ^- 


11 


-2^ - 


^ B 




^2B" 




33 — 


-R22 sin^ 


6* 



A 



32A'2 
' A4 



,A' 



,B^'^ 
'B2 



I 9 A^ I nB^ _ B'A' 



BA 

, 9„2 B'2 , 2 B'A' 



dB 
dr 



(4.6) 



(4.7) 
(4.8) 

(4.9) 



where A' = f, 

Second Order Skew Tensors 

Using the tetrad (4.2), (4.3), the metric(4.4), (4.5) and using the definition given 
in Table(l), we found that all components of the 2nd order skew tensor Ea/s vanish iden- 
tically, as expected. Consequently, the skew part of the field equations (3.9) is satisfied 
identically. 
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Scalars 

The scalars which are necessary for the field equations (3.3) can be obtained from 
the above components of the second order symmetric tensors. These are, 

^ = ^'(7? + 2^J (4-10), 

and 

-/I o ^4 o -i4 y4. n ,o // 8 / , ^ 

R = 2 ^ 25^— + 4B^-nr - 4— 5^ + 65'^ - 4S"5 - -B'B. (4.11 

The Field Equations 

Using the above calculated tensors and the field equations (3.1) in its mixed form, we 
get the following set of differentional equations, 

S'2 B' A'B' A" A'^ A' ^ 

B'^ B' B'A' A' 

„2r B" 5'2 B' A'B' A" A'^ A' ^ 

Assuming that B is non- vanishing, to prevent singularities in (4.3) and (4.4), then we can 
write the above set of differentional equations in the form, 

'2 M 

(4.12) 

(4.13) 
(4.14) 



4^- 

B 


B' 

^B^^^B-r- 


A'B' 


-4 




A' 

Ar 


= 




B'^ B' 


,B'A' 

-\a 


A^ 


A' 


= 0, 




4- 


B'^ B' 
4— + 4 




-4 


A' 


A 
Ar 


= 


From (4.12) and (4.14) 


we get 




B' 









(4.15) 

which gives, by integration 

^-^.^ (4-16) 

where C* is the constant of integration. The solution (4.16) satisfies equation (4.13) 
without any further condition. So, we see that the set of differentional equations (4.12)- 
(4.14) is not sufficient to determine the explicit forms of the two function A and B. This 
will be discussed in section 5. 

In what follows we are going to obtain some physical information from the geometric 
model given above. 
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4.1 A Method for Fixing the Unknown Function 

In this theory, it is clear that the material-energy tensor is a pure geometric object (3.6), 
which has the definition 

^* def 1 , 1 , 1 

Using the second order symmetric tensors calculated above, we can get the following 
non- vanishing components of the geometric material-energy tensor (3.6), 



t; 



00 



l2 



T* — r'^ 
-'-22 — ' 



'A'B' 


A" 


3A'^ 
^ 2 






AB 


~ ~A 


Ar 




T*,= 




B' 


A'^ 
2^2' 




B" 


+ 


B' 
-B-r + 


A'B' 


A'^' 




AB 


2A2 



and 



r3*3 = t;^ sin^ e. 



(4.17) 

(4.18) 
(4.19) 



Using the definition T*^^ =^ g'^^T*^ ,(4-15) and the metric (4.5) we can write the above 
components in the form. 



rj-\*0 j^2 



B" B' 

-3— + 2 h4 — 

52 B Br 



5'2 B' 

3^ - 2 — 
52 



B" B' 

~b' " 'b^~^^ 



(4.20) 

(4.21) 
(4.22) 



and 



. 3 ~ ^ .2 



If we assume that we have a perfect fiuid, then the geometric energy-momentum tensor 
will satisfy the following relation 

T*\ = = r*% . (4.23) 

Then, using this relation, i.e. equating (4.21) and (4.22), we get 

B" B' 
-B^^n^-Wr^''^ 

which is a second order differentional equation in [B) only. Integrating this differentional 
equation twice we get 

^2 

(4.24) 



B = {IC- + 7Ci) 



So, from (4.16) we can write 



- ^ . . (4.25) 



5' (7Cf + 7Ci) 
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This shows that the field equations (3.3) can fix the unknown functions of the model only 
if we attribute some properties to the material distribution. 



4.2 The Metric of the Associated Riemannian Space 

In order to get some physical information about the solution obtained, (4.24) and (4.25), 
we are going to write the metric of Riemannian space, associated with the AP-space (4.3), 
which is given, in general, by 

ds^ = g^^dx^'dxr (4.26) 
Substituting the solution obtained, (4.24) and (4.25), into (4.4) we can write, 

+ ^ Adr^ + r^de^ + r2 sin^ Qd^'') ' ^^•2'^) 

in order to facilitate comparison with the results of GR, we are going to write (4.27) 
assuming that (C*)^ — —C2 ( where is another arbitrary constant) and dr — ids, then 
we get, 

dr2= ^(7c4 + 7Ci)ldt2 

^ M'' + r^de^ + r^sm^9dip^) ' (^"^S) 

(7C^+7Ci)7 ^ ^ ^ 

which is the metric of the pseudo-Riemannian space associated with the AP-space, with 
the same signature used usually in such applications of GR. 

Boundary Conditions 

The non-vanishing of some of the components of the geometric material-energy 
tensor (3.6), shows clearly that the solution obtained is an interior one, having spherical 
symmetry. To fix the constants of integration C, Ci and C2, we are going to make a 
matching between the metric (4.28) and the Schwartzchild exterior metric [2] at the 
boundary r = a. Now, let us define the geometric density (p*o) and the geometric pressure 
(p*o) (these quantities are measured in relativistic units), we get 

pMr)"^ -T*\ , pVr) T*^ = = T*%. (4.29) 

Then we can determine from equation (4. 21), (4. 29) using the solution (4.24) and 

(4.25). This can be written in the form. 
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Then, assuming that p^r) = at the boundary r = a we get 

7Ci 



2a2 ' 



and from matching the metric (4.28) with the Schwartzchild exterior solution, which can 
be written in the form [2] 

at r = a we get, 

4 1 

7Ci = 



3^(1 + 1)^^ ' 
and 

c. 



Consequently we can write the solution (4.24) and (4.25) in the form. 



B = ^^77^^^ , (4.31) 

and 



Then the metric (4.28) can be written in the form, 

-2 _ I, /I 7 „2n| 7.2 ^2 ,,J2 , 2in2 , „2 •„2 



(4.32) 



dr' = 6i(l - -^r')^(ir -4f^{dr' + r^d^ + r'sin' Odip') , (4.33) 



where 



, def /~4 4 (1 — ^)^ 

and 

62 (.ii) 

Then substitute from (4.29), (4.31) and (4.32) into (4.20) and (4.21) using (m) and {vii) 
we can evaluate PQ{r) and PQ{r) in the form, 

Pl^r) = 4r J^~.^^L^ > (4-34) 



a 



and 



^S(r)=3-^-^^-^. (4.35) 

These two equations give the values of the geometric pressure and density as functions of 
the radial distance r. 
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5 Discussion 



In the present work, we have apphed a pure geometric field theory, in which the tensor 
representing the material-energy distribution, given by (3.6), is a part of the geometric 
structure used. The field equations of the theory are applied to a geometric structure 
with spherical symmetry. It appears that the field equations are not sufficient to fix the 
unknown functions in this case. The non- vanishing of the components of the geometric 
material-energy tensor, as clear from equations (4.17), (4.18) and (4.19), shows that the 
model represents the field within a material distribution. We show that, the unknown 
functions of the model can be fixed by assuming, in addition to the field equations, that 
a perfect fiuid is filling the spherical configuration. The geometric density and pressure 
are defined (4.29). Some boundary conditions are used to fix the constants of integration. 
This is done by assuming, at the boundary {r — a), that: 

(i) The geometric pressure (pq) vanishes. 

(ii) The metric (4.28) obtained from the solution and the Schwartzchild 
exterior solution should be the same at this boundary. 

Now, we can direct the attention to the following points: 

(1) The field theory used can be classified as a pure geometric one. In other words, 
all physical quantities are represented by geometric objects that are constructed from the 
building blocks of the geometric structure used, the AP-structure of the Riemannian- 
Cartan type. This is done to avoid the use of a phenomenological material-energy tensor 
and an equation of state. 

(2) The geometric structure with spherical symmetry, given in section 3, used in the 
present work, is that usually used in such applications (cf.[10],[ll]). The non- vanishing of 

2nd-ordcr skew tensors (see Table 1) may give rise to interactions other than gravity. So, 
to restrict the model to gravitational interactions only, we have assumed the vanishing of 
the function D which leads to the vanishing of 2nd order skew tensors, relevant to the 
theory used. 

(3) It is shown that the field equations (3.1) are not sufficient to determine the un- 
known functions A{r)^ B{r) of the model. Instead, it gives a relation (4.16) between A 
and B. Any pair of functions satisfying this relation represents a solution. Then the 
relation (4.16) implies a family of solutions. On the other hand, the spherically sym- 
metric geometric structure, used in the present study, can be considered as admitting 
representation of several physical situations. Once one of such situations (perfect fluid) is 
assumed, the functions can be fixed as shown in (4.31), (4.32). This has to be expected 
from the beginning, since the components of the material-energy tensor (4.20), (4.21) and 
(4.22), are pure geometric objects, without imposing any restrictions, expect for spherical 
symmetry. 

(4) The boundary condition used in the previous section, r = a, can be considered as 
a boundary between two regions, the first is characterized by 7^ while the second is 



12 



characterized by — 0. We can consider the value (a) as representing the radius of the 
spherical configuration ( a static or a slowly rotating star). 

The radial distance (r), measured from the center of the star, can be written in terms 
of the dimensionless quantity, 

def r 

q ^ - ■ 
a 

In this case the solution (4.31) and (4.32), the geometric pressure (4.34) and the geometric 
density (4.35) can be written in terms of q and the parameters bi and 62, of the previous 
section, as 



B^^-^^^, (5.1) 



-'2 



I 



^ = ^ 2 , (5.2) 



a 



and 



„2 , 



1 fl-^) 

Po(r)=3— --^-y^. (5.4) 
a^02 (1 - jq^)~ 



(5) Corona: If we evaluate the geometric density (4.35) at r = a, or g = 1, we get 

, 0.7771717 
Po=^^ (5.5). 

Imposing the condition for vanishing density (pq = 0) on (5.4), we get from (4.35) 

q^V2. (5.6) 

So, we have a region (spherical shell), outside the star, within which Pq — while ^ 0. 
The width of this region is in the range -\/2 > g > 1. This region can be considered 
as a corona of the star. However, the corona of the Sun is not an exact spherical shell. 
This is due to the magnetic activity of the sun. Recalling that the model in the present 
treatment is a pure gravity model, i.e. without considering any electromagnetic influence. 
This may give an interpretation for the difference between the shape of the corona, given 
in the present work, and that of an actual star. 

(6) Equation of state: Using the definitions of the geometric pressure (5.3) and density 
(5.4), we can write the following relation between (pg) and (pq), as 

pI = ii^^ pI . (5.7) 
° 3(1 -4) ° ^ ' 

This represents a general equation of state that is obtained from ,and not imposed on, 
the model. Equation (5.7) is one of the advantages of using a pure geometric theory 
of gravity. In fact equation (5.7) represents a continues spectra of characteristics of the 
material distribution within the star. This spectra can not be obtained if one uses the 
conventional approach, since the equation of state is imposed on the model, in that case. 
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The relation between the geometric material-energy tensor and the phenomenological 
one (written in cgs unit) is obtained by comparing the R.H.S. of (3.5) and the R.H.S. of 
(1.1) so, we can write 

T*^^ = -K T^. (5.8) 

Then we get, 

* ' 
Po Po 

where po and po are the physical pressure and density of the material distribution, re- 
spectively. The above relation and Figure 1 show how {^) varies from the center to the 
surface of the star. 
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(7) Radiation Zone: It is clear from the equation of state (5.7) (or from Fig.l) that as 
q 0. — ^ ^. This value characterizes the situation of a radiation dominant central 
zone. This situation washes out gradually as we move towards the surface of the star 
(g = 1) , at which pq = 0, which characterizes energy transfer in a class of main sequence 
stars [12], according to Schwarzschild classification. 

(8) Singuleirities and Convection Zone: It is clear that the solution given (5.2) is 
singular at g = ~ 0.755928. It is clear from that (5.3)) and (5.4) that the geometric 
pressure and density are also singular at the same value of q , while the equation of state 
(5.7) is regular. 

The appearance of singularity, may indicate that the model is applied outside its 
domain of applicability. In other words, the model may include one or more assumptions 
that is not applicable at the singularity. Let us recall the case of the Sun. There is a 
convection zone at a depth about 28% of the radius of the Sun, i.e. 72% of the radius 
as measured from the center of the Sun [13]. The perfect fluid assumption within the 
convection zone is no longer valid. The presence of singularity in the present work may 
indicate the existence of a zone in which perfect fluid assumption is violated, a convective 
zone. Fig. 2 shows the variation of the density and pressure as function of q. 
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Figure 2: The relation between q and Pq is represented by the dashed line, while the 
relation between q and p*Q is represented by the solid line 



(9) The model obtained, in the present work, is far from representing a compact 
object. This is because the equation of state for such objects is the adiabatic one [14], 
given by 

po = Kpl^- , (5.9) 

where {K) is constant and (n) is the polytropic index. This index varies from (|) for 
non-relativistic particles to (3) for relativistic particles. It is clear that the equation of 
state obtained (5.7) has a polytropic index tending to infinity, which characterizes an 
isothermal configuration [12]. 

6 Concluding Remarks 

The situation here is similar to the Schwartzchild exterior solution in which pure geometric 
considerations led to many successful physical interpretation and predictions. This shows 
that Einstein geometerization scheme can be extended to cover domains with material 
distribution. 
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We may conclude that the model obtained in the present work represents a lower main 
sequence star with a corona, radiation zone and convection zone. 

Finally, we would like to point out that the introduction of electromagnetic effects, on 
the model, may throw more light on the physical situation within the star. This can be 
done by relaxing the condition, D{r) = 0, imposed on the model. 
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